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In the framework of the theory of open systems based on completely positive quan- 
tum dynamical semigroups, we solve in the asymptotic long-time regime the master 
• equation for two independent harmonic oscillators interacting with an environment. 

We give a description of the continuous-variable asymptotic entanglement in terms 
of the covariance matrix of the considered subsystem for an arbitrary Gaussian input 
state. Using Peres-Simon necessary and sufficient condition for separability of two- 
mode Gaussian states, we show that for certain classes of environments the initial 
state evolves asymptotically to an entangled equilibrium bipartite state, while for 
^ ' other values of the coefficients describing the environment, the asymptotic state is 

^ ■ separable. We calculate also the logarithmic negativity characterizing the degree of 

O^, entanglement of the asymptotic state. 

> : 1 Introduction 
in 
i> 

■ When two systems are immersed in an environment, then, besides and at the same time 
^ ! with the quantum decoherence phenomenon, the environment can also generate a quantum 
^ I entanglement of the two systems [HIS]. In certain circumstances, the environment enhances 
I the entanglement and in others it suppresses the entanglement and the state describing the 
O ■ two systems becomes separable. The structure of the environment may be such that not 
only the two systems become entangled, but also such that the entanglement is maintained 
^ , for a definite time or a certain amount of entanglement survives in the asymptotic long-time 
regime. 

In the present paper we investigate, in the framework of the theory of open systems 
based on completely positive quantum dynamical semigroups, the existence of the continu- 
ous variable asymptotic entanglement for a subsystem composed of two identical harmonic 
oscillators interacting with an environment. We are interested in discussing the correlation 
effect of the environment, therefore we assume that the two systems are independent, i.e. 
they do not interact directly. The initial state of the subsystem is taken of Gaussian form 
and the evolution under the quantum dynamical semigroup assures the preservation in time 
of the Gaussian form of the state. We only investigate here the asymptotic behaviour of 
the subsystem states. The time evolution of the entanglement, in particular the possibility 
of the so-called "entanglement sudden death", that is suppression of the entanglement at 
a certain finite moment of time, will be discussed in a future work. 

The organizing of the paper is as follows. In Sec. 2 we write the equations of motion in 
the Heisenberg picture for two independent harmonic oscillators interacting with a general 
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environment. With these equations we derive in Sec. 3 the asymptotic values of variances 
and covariances of coordinates and momenta which enter the asymptotic covariance matrix. 
Then, by using the Peres-Simon necessary and sufficient condition for separabihty of two- 
mode Gaussian states [3i |4j, we investigate the behaviour of the environment induced 
entanglement in the limit of long times. We show that for certain classes of environments 
the initial state evolves asymptotically to an equilibrium state which is entangled, while for 
other values of the parameters describing the environment, the entanglement is suppressed 
and the asymptotic state is separable. The existence of the quantum correlations between 
the two systems in the asymptotic long-time regime is the result of the competition between 
entanglement and quantum decoherence. A summary is given in Sec. 4. 



2 Equations of motion for two independent harmonic 
oscillators 

We study the dynamics of the subsystem composed of two identical non-interacting os- 
cillators in weak interaction with an environment. In the axiomatic formalism based on 
completely positive quantum dynamical semigroups, the irreversible time evolution of an 
open system is described by the following general quantum Markovian master equation for 
an operator A (Heisenberg representation) [51E]: 

^ = 'j-[H,A{t)] + ^E(^'K^)'^.] + [vl^mm. (1) 

Here, H denotes the Hamiltonian of the open system and Vj, Vj which are operators 
defined on the Hilbert space of H, represent the interaction of the open system with the 
environment. Being interested in the set of Gaussian states, we introduce those quantum 
dynamical semigroups that preserve that set. Therefore H is taken to be a polynomial of 
second degree in the coordinates x,y and momenta Px,Py of the two quantum oscillators 
and Vj,Vj are taken polynomials of only first degree in these canonical observables. Then 
in the linear space spanned by the coordinates and momenta there exist only four linearly 
independent operators V^=i, 2,3,4 [Z] : 

Vj = a^jPx + ayjVy + bxjX + byjy, (2) 

where axj,ayj,bxj,byj G C. The Hamiltonian H of the two uncoupled identical harmonic 
oscillators of mass m and frequency u is given by 

H = ^ipl+pl) + ^i^' + y^ (3) 

The fact that the evolution is given by a dynamical semigroup implies the positivity of 
the following matrix formed by the scalar products of the four vectors a^, b^, a^, by whose 
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entries are the components axj,bxj,ayj,byj, respectively: 



-n 



^ ^x^x) i^x^x) i^x^y) i^x^y) \ 

(Ka-x) (KK) (Kay) {hxhy] 

(a^a^) (ayhx) i^y^y) i^y^y] 

V (bj^a^.) (bj^b 

x) O^y^y ) (WW / 



(4) 



For simplicity we take this matrix of the following form, where all coefficients D^x, Dxp^,.. 
and A are real quantities: 



V 



Dxx 

-Dxp, + ihX/2 



-Dxp^ - ihX/2 

Dr 



'xy 



'PxPx 

-a 



-D 



XPy 



n 



ypx 

PxPy 



D 
-R 
D 



xy 

ypx 



yy 



-D 



ypy 



ihX/2 



Dxpy \ 

^PxPy 

-Dypy - ifiX/2 
D 



(5) 



'PyPy 



It follows that the principal minors of this matrix are positive or zero. From the Cauchy- 
Schwarz inequality the following relations for the coefficients defined in Eq. hold (from 
now on we put, for simplicity, ^ = 1): 



D D — > Q D D — 

^xx-^iiii xy — ' xx-^pxPx , 

D D - D'^ > D D - D"^ 



> ^, DxxBp^p., - Z^L, > 0, 
A2 



> 



Dt 



yy-^PxPx ^ypx — ^yy^PvPy yPy — ^ ' PxPx^PyPy ^PxPy 



> 0. 



(6) 



The covariance of operators Ai and A2 can be written with the density operator p, 
describing the initial state of the quantum system, as follows: 



^A,A,{t) = -Tt{p{A,A2 + A2A,){t)). 
We introduce the following 4x4 bimodal covariance matrix: 



/ CTx 



a{t) 



XPx 



a. 



xy 



XPx 



ex. 



PxPx 



(J, 



ypx 



0"; 

a. 



XPy 



a. 



xy 



a. 



ypx 



yy 



Xpy ^PxPy ^ypy 



PxPy 

^ypy 

^PyPy I 



By direct calculation we obtain [7] (T denotes the transposed matrix): 

do 



dt 



Ya + aY' + 2D, 



where 



Y 
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-A 


1/m 









^ Dxx 


Dxpx 


Dxy 


Dxpy \ 




—muj'^ 


-A 








, D = 


Dxpx 


DpxPx 


Dijpx 


DpxPy 










-A 


1/m 




Dxy 


^ypx 


^yy 


^ypy 


V 








— mcj^ 


-A j 




\ Dxpy 


DpxPy 


^ypy 


^PvPy 1 



(7) 



(8) 



(9) 
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The time-dependent solution of Eq. ([9]) is given by [7] 

a{t) = M(t)(a(0) - a(oo))M^(t) + a{oo), (11) 

where M{t) = exp{Yt). The matrix M{t) has to fulfil the condition lim^^oo (t) = 0. In 
order that this limit exists, Y must only have eigenvalues with negative real parts. The 
values at infinity are obtained from the equation [7] 

Ya{oo) + a{oo)Y^ = -2D. (12) 



3 Environment induced entanglement 

The two-mode Gaussian state is entirely specified by its covariance matrix a (IE]), which is 
a real, symmetric and positive matrix with the following block structure: 

a=(^T (13) 

where A, B and C are 2x2 matrices. Their entries are correlations of the canonical 
operators x,y,Px,Py, A and B denote the symmetric covariance matrices for the individ- 
ual reduced one-mode states, while the matrix C contains the cross- correlations between 
modes. In the following we shall consider environments for which Dxx = Dyy,Dxp^ = 
Dypy, Dp^p^ = Dpypy, Dxpy = Dyp^. Then both unimodal covariance matrices are equal, 
A = B and the entanglement matrix C is symmetric. With the chosen coefficients, we 
obtain from Eq. (fT2l) the following elements of the asymptotic entanglement matrix C: 



( ^ /X -m'^uj'^Dxy + 2m\Dxpy + Dp^p 
axp^ioo) = aypAoo) = 2m(A^^a;^) ' ^'^^ 



, , rrruj'^Dxv — 2miu^XDxv + (2A -|- lu )Dr, r, , ^ 

(J fool = — ■ - ^^'^ flfil 



The elements of matrices A and B are obtained by taking x = ?/ in the previous expressions. 
We calculate the determinant of the entanglement matrix and obtain: 

Gaussian states with det C > are separable states, but for det C < it may be possible 
that the asymptotic equilibrium states are entangled. In order to investigate whether the 
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environment can actually entangle the two independent systems, we can use the partial 
transposition criterion [21 H]: a state is entangled if and only if the operation of partial 
transposition does not preserve its positivity. Simon |4| obtained the following necessary 
and sufficient criterion for separability: S > 0, where 

S = det AdetS+ (| -\detC\f - Ti[AJCJBJC"^ J] - ^(det A + det 5) (18) 

and J is the 2x2 symplectic matrix. 

In order to analyze the possible persistence of the environment induced entanglement 
in the asymptotic long-time regime, we consider the environment characterized by the 
following values of its parameters: rn?uj'^Dxx = Dp^p^, D^p^ = 0, m^ui'^Dxy = Dp^p^. In 
this case the Simon expression f|T8l) takes the form: 

S={ ^^f ^ + - 7 - 4^^7TT^^^. 19 

\ A2 A2 + cj2 4/ A2(A2 + u;2) 



For environments characterized by such coefficients that the expression ([T9|) is strictly 
negative, the asymptotic final state is entangled. In particular, if D^y = 0, we obtain that 
S < 0, i.e. the asymptotic final state is entangled, for the following range of values of the 
coefficient D^py characterizing the environment ^J: 

A 2<v/pT^< A +2' ^^^^ 

where the coefficient D^.^^ satisfies the condition mwDxx/^ > 1/2, equivalent with the 
unimodal uncertainty relation. If the coefficients do not fulfil the inequalities (1201) . then 
5* > and the asymptotic state of the considered system is separable. 

We apply the measure of entanglement based on negative eigenvalues of the partial 
transpose of the subsystem density matrix. For a Gaussian density operator, the neg- 
ativity is completely defined by the symplectic spectrum of the partial transpose of the 
covariance matrix. The logarithmic negativity E = — | log2[4/((T)] determines the strength 
of entanglement for E > 0. U E < 0, then the state is separable. Here 



/(a) = ^(det A + det B) - det C 



1 ^2 



2 (det A + det B) - det C 



det a. (21) 



In our case the logarithmic negativity has the form 



E = - log2 



A VA^ 



(22) 



It depends only on the diffusion and dissipation coefficients characterizing the environment 
and does not depend on the initial Gaussian state. In Fig. 1 logarithmic negativity E is 
represented versus scaled diffusion coefficients D^x and D^py- 
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Figure 1: Logarithmic negativity versus scaled diffusion coefficients. 

4 Summary 

In the framework of the theory of open quantum systems based on completely positive 
quantum dynamical semigroups, we investigated the existence of the asymptotic quantum 
entanglement for a subsystem composed of two uncoupled identical harmonic oscillators 
interacting with an environment. By using the Peres-Simon necessary and sufficient con- 
dition for separability of two-mode Gaussian states, we have shown that for certain classes 
of environments the initial state evolves asymptotically to an equilibrium state which is 
entangled, i.e. there exist non-local quantum correlations for the bipartite states of the 
two harmonic oscillator subsystem, while for other values of the coefficients describing 
the environment, the asymptotic state is separable. We determined also the logarith- 
mic negativity characterizing the degree of entanglement of the asymptotic state. Due to 
the increased interest manifested towards the continuous variables approach to quantum 
information theory, these results, in particular the possibility of maintaining a bipartite 
entanglement in a diffusive-dissipative environment for asymptotic long times, might be 
useful for applications in the field of quantum information processing and communication. 
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